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ABSTRACT 


Wireless communication is employed to connect mobile computers in a networked 
environment for information exchange. In a tactical space, sensors and computers 
typically need to operate on-the-move while transmitting data over both short and long 
distances in different terrain and conditions. The wireless communication is thus 
susceptible to effects of Doppler shift and channel fading. In addition, when security and 
anti jamming features are required, such as frequency-hopping techniques, then coherent 
signal detection is difficult and noncoherent modulation is used instead. 

Our study will focus on the bit error rate (BER) performance analysis of 
noncoherent orthogonal modulation, specifically M-ary frequency-shift keying (MFSK) 
and code-shift keying (CSK) modulation, in both additive white Gaussian noise (AWGN) 
and for a Rayleigh fading channel with Doppler shift. The potential applications include 
communications between mobile computer-sensor devices, such as a mobile ground 
control station maintaining a datalink with UAV. 
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EXECUTIVE SUMMARY 


Wireless communication is employed to connect mobile computers in a networked 
environment for information exchange and is susceptible to effects of Doppler shift and 
channel fading. Examples of potential applications include ground station datalink with 
UAV and satellites relaying information to aircraft. In this thesis, bit error rate (BER) 
performance analysis for noncoherent demodulation of M-ary frequency-shift keying 
(MFSK) and code-shift keying (CSK) in both Rayleigh fading channel and additive white 
Gaussian noise (AWGN) channel with Doppler shift was carried out. 

Numerical analysis was done using Matlab followed by simulation of MFSK with 
Simulink. In general, MFSK performed better than CSK by approximately 2 dB to 

4.4 dB at BER of about 10 4 for a Rayleigh fading channel and about 1.8 dB to 3.4 dB 
for an AWGN channel when considering a Doppler shift f d normalized to the symbol 
duration T s at f d T s = 0.3. This was due to Walsh coded signals in CSK becoming 

unrecognizable at the demodulator when they experienced Doppler shift, especially for 
higher modulation orders. 

It is possible to improve BER performance for both MFSK and CSK by 
increasing the order of modulation; however, there are diminishing returns and it is not 
recommended to go higher than M = 8 since continued improvement is marginal. A 
much higher SNR is required for a Rayleigh fading channel to achieve the same BER as 
for no fading. The case was presented for BER of 10 4 and Doppler shift of f tj T = 0.3 . 
For CSK the signal-to-noise ratio (SNR) difference was approximately 26 dB to 29 dB , 
while MFSK required a SNR of about 17 dB to 28 dB higher for a Rayleigh fading 
channel. 
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I. 


INTRODUCTION 


A. OVERVIEW OF WIRELESS COMMUNICATION 

Wireless communication systems have ubiquitous presence in both military and 
commercial applications due to significant operational advantages for long range 
operations where wires cannot be run and easily maintained. Established industry 
standards, such as WiFi (IEEE 802.11), permit the setup of communications quickly 
when existing infrastructure is not readily available or accessible, for example in areas of 
conflict zones or humanitarian assistance and disaster relief (HADR) missions where 
foreign aid can be rendered at a moment's notice. 

Wireless communication also reduces the cost of network deployment and 
expansion. Hardware costs have dropped precipitously, making computers widely 
available to support the growth of digital content on the Internet. Wireless connectivity 
has increased with long-term evolution (LTE) and worldwide interoperability for 
microwave access (WiMAX) as the emerging standards for delivery of last mile wireless 
broadband access, potentially enabling mobile Internet access worldwide. 

For the military, a similar trend is evolving. The concept of network centric 
warfare is key to future combat systems (FCS) development and will allow combatants, 
logisticians, weapon and sensor systems to upload, download, and share real-time 
pertinent information about their mission. Examples include soldiers tapping into C4I 
networks with mobile computers for situational updates on the battlefield, logisticians 
downloading mission support, and operators maintaining datalinks between an unmanned 
aerial vehicle (UAV) conducting surveillance missions, the ground control station (GCS), 
and orbiting satellites (Figure 1). 
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Figure 1. A typical UAV operating scenario where datalinks are established between 
the ground control station (GCS) and satellite in orbit (From [1]). 

Mobile computing will continue to grow in usage, and wireless communication 
will form the backbone of both the Internet and military C4I networks (Command, 
Control, Communications, Computers, and Intelligence). 

B. LIMITATIONS TO PERFORMANCE BY FADING AND DOPPLER 

The performance of wireless communication is limited by the environment 
because the transmitted signal must propagate across the mobile radio channel to arrive at 
the receiver. In addition to the background noise and interference from other sources in 
the medium the signal must contend with, there may be no direct line-of-sight (LOS) path 
between the transmitter and receiver in typical urban areas and terrain with mountains 
and forests. The signal must rely on non-free-space propagation modes: reflection, 
diffraction and scattering [2]. 

The received signal is a combination of many copies of the transmitted signal that 

have propagated through multiple paths in the medium, with each path having its own 

attenuation, phase shift, and time delay on the signal. Arriving copies of the transmitted 
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signal at the receiver may interfere either constructively or destructively with one 
another. Overall, the received signal becomes attenuated by multipath fading. 

Mobile computers can also communicate wirelessly while on the move, so there is 
some relative motion between transmitter and receiver, such as the UAV with the GCS 
and orbiting satellite (Figure 1). Even for a given set of coordinates for the transmitter 
and receiver, non-static objects in the medium can function as moving reflectors. This 
causes the received signal's frequency to change and results in the loss of orthogonality, 
an important property that uniquely defines symbols. The effect of Doppler shift can 
severely degrade the performance of orthogonal modulation schemes such as M-ary 
frequency-shift keying (MFSK) and code-shift keying (CSK). 

C. THESIS OBJECTIVES 

The performance of noncoherent orthogonal MFSK in the presence of timing and 
frequency errors is assessed in additive white Gaussian noise (AWGN) for modulation 
orders of 2, 4, and 8 in [3]. This is with respect to bit error performance due to the lack 
of perfect time and frequency synchronization, and the effect of fading is not considered. 
In addition, CSK performance with interleaving is evaluated in fading channels for 
modulation orders of 4, 8, 16, and 32 in [4] by approximation using the Chernoff bound. 
In these references, however, the effect of Doppler shift is not included. 

In this thesis, the bit error performances of MFSK and CSK in the presence of 
Rayleigh fading channel and Doppler shift will be studied. The performance of MFSK 
and CSK in AWGN with Doppler shift will also be included. Using Matlab, numerical 
analysis will be carried out to evaluate the bit error rates (BER) with respect to the bit 
energy-to-noise ratio E b / N o . In general, a BER of no greater than 10 1 to 10 5 (about 
one bit error per 10,000 bits) is required for reliable digital data communication; 
therefore, the selection of E b / N a for the BER performance charts of MFSK and CSK 
will cover these ranges of BER. Simulation using models developed in Simulink will 
also be carried out and compared to results from the numerical analysis. 
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D. THESIS OUTLINE 


This thesis is organized into seven chapters. The introduction to the thesis is 
given in Chapter I. The background to Doppler effect and the Rayleigh fading channel is 
covered in Chapter II. The numerical analysis of noncoherent demodulation of MFSK 
and CSK for a Rayleigh fading channel with AWGN and Doppler shift are presented in 
Chapters III and IV, respectively. In Chapter V, the discussion is extended to include 
BER performance of MFSK and CSK in AWGN with Doppler shift. The simulation 
results from Simulink are presented and compared to numerical analysis in Chapter VI. 
Finally, the conclusions based on results in Chapters III through VI are stated in Chapter 
VII. 
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II. BACKGROUND 


A. DOPPLER EFFECT 

In wireless communications, radio waves of different frequencies are employed to 
transmit energy from a transmitter to an intended receiver. Doppler effect (or Doppler 
shift) occurs when there is relative motion between the transmitter and receiver. The 
perceived wavelength of the transmitted radio wave by the receiver will appear to change, 
causing the received frequency to differ from the original sent frequency. While sensor 
systems, such as Radar with pulse compression technology, may exploit Doppler shift to 
determine speeds of moving targets, the Doppler effect has far more undesirable 
consequences for typical wireless communications systems. 

When the receiver is moving towards the transmitter with a radial velocity ofv, 
the Doppler shift f d is 

f v 

f d =^~ ( 1 ) 

c 

where f c is the carrier frequency and c is the speed of light. The value of f d is negative 
if the receiver is moving away from the transmitter. Instead of detecting the carrier 
frequency f c at the receiver, the received frequency / will be 

fr=fc+f d ={l + ~)fc- ( 2 ) 

^ c) 

As a result, the received frequency is higher than the transmitted frequency. 
Conversely, the frequency is lower if the receiver is moving away from the transmitter. 
Because the frequency is altered, the received signal loses orthogonality, which is 
important for its detection of noncoherent schemes such as MFSK and CSK. With 
Doppler shift, receiver performance with orthogonal signaling degrades (i.e., the BER 
increases). 
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B. THE RAYLEIGH MULTIPATH FADING CHANNEL 


The mobile radio channel is a multipath fading channel that causes attenuation, 
phase shifting, and time delays in the received signals. The channel exhibits a filtering 
nature and may be modeled as a linear filter with a time varying impulse response. The 
time variation is due to relative motion between the receiver and transmitter and 
reflectors in the medium. Similar to the filter taps in a finite impulse response (FIR) 
filter, every path in the multipath radio channel can be represented as a channel tap. 
Each of these paths is called a resolvable path , which has an arbitrary attenuation, phase 
shift and path delay on the transmitted signal. 

The transmitted bandpass signal s(t) is represented in the complex envelope 
representation as 

s(t) = Re{s L (t)e j2 * f '‘} (3) 

where s L (t ) is called the complex envelope and the carrier frequency is /' . As the 
transmitted signal s(t) is propagated through different resolvable paths in the fading 
channel, multipath copies of the signal are generated. 

Assuming there are i different resolvable paths, we see that each corresponding 
multipath copy is itself the sum of all the uncorrelated signals having similar (or 
approximately the same) path delay r j (t) but having k different attenuations a ik (t) and 

phase shifts <j> ik (t ). The Doppler shift of the path is also included in (fi k ( t ). 


The paths that uncorrelated signals take are called non-resolvable paths and are 
grouped together in a single multipath copy since they have the same path delay zv(f). 


Since the channel’s filtering nature is time-varying, these multipath channel parameters 
are functions of time. The received signal is now represented as 


x(t) = Re J X Z a ik 0 ) s l ( t ~ T , (t)) e 


j[2nf c (t-T,(t))+0j k (t)] 


i k 


= Re< 


X2 l a ik (t)s L {t-^(t))e 




A 


(4) 


V i k 


yfinfet 
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The complex envelope x L (t) is then expressed as 


x l (0 = ZE%(0e" j[2 " /<ri( ' ) ' 4( ' )1 ^ (r —r # (r)) 

i k 

= Y J h i(t)s L ( t ~T i (t)) 


(5) 


where /i ; (?) is called the i"‘ complex path attenuation or the i' h channel tap. Since there 
are many multipaths in the channel, the central limit theorem can be applied to model the 
channel tap lift) as a complex Gaussian random variable [5]. 

If a dominant line of sight (LOS) path exists between the transmitter and receiver, 
the magnitude \hft)\ of the i !h channel tap follows a Rician distribution. In the presence 

of many obstacles, such as in a typical dense urban or forested environment, the 
multipath scatterers are relatively independent with similar probability distribution, and 
\hft)\ will then follow a Rayleigh distribution. 


C. CHAPTER SUMMARY 

In this chapter, it was observed that Doppler shift f d caused by relative motion 

between the transmitter and receiver alters the perceived frequency of incoming signals at 
the receiver. This affects the orthogonality of the signals and is detrimental to the bit 
error rate (BER) performance of orthogonal signaling such as MFSK and CSK. A 
multipath fading channel causes the transmitted signal to experience attenuations, phase 
shifting, and time delays. Multiple copies of the signal interfere with one another at the 
receiver. In this thesis, the Rayleigh fading channel is considered, and its filtering nature 
is modeled with the channel tap \h k |. 

In the following chapters, the performance of noncoherently detected MFSK and 
CSK are evaluated with the effect of Doppler shift. In particular, consider the case where 
the channel is assumed to have flat fading with one channel tap \h k \ and is slow fading. 
The first assumption implies that all frequency components of the signal experience the 
same magnitude of fading, \h k \. The second assumption implies that the channel phase 

shift 6 k is considered constant within the symbol duration T s . 
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III. PERFORMANCE OF M-ARY FREQUENCY SHIFT KEYING 
WITH NONCOHERENT DETECTION 


MFSK is an orthogonal signaling scheme which uses a set of M = 2 k sinusoids to 
represent M distinct symbols, each of k bits. For example, a 64-FSK scheme uses 64 
orthogonal sinusoids for 64 different six-bit symbols. Two real-valued sinusoids s, (t) 

and s 2 (t) are orthogonal to each other over the interval [a,b] when their inner product is 
zero as in 

rb 

sdt)s 2 (t)dt = 0. (6) 

J a 

At the demodulator, the received signal is typically passed through a bank of M 
matched filters with each branch containing one of the M number of transmitted 
waveforms generated locally as reference signals. For coherent demodulation , these 
reference signals must be matched exactly in frequency and phase to the received signal 
using techniques such as phase lock loops. The received signal is identified at the branch 
where the output of the integrator is maximum. For example, if s^t) was received, the 

peak value will occur only at the integrator output in the corresponding branch with v, (t) 
as the sample waveform to be detected. From Equation (6), all the other branches should 
yield a relatively low output since their sample signals s 2 (t),...,s M (t ) are orthogonal to 

s x (t) with inner products close to zero. Thus, maintaining signal orthogonality is 
important for correct signal detection. 

Cases exist where the phase of the received signal cannot be recovered accurately 
or tracked easily for coherent demodulation, such as in a fast fading Rayleigh channel or 
a communication system using frequency-hopping spread spectrum. Additionally, 
coherent detection does not work well if the transmitted signal is phase-shifted by the 
radio channel during propagation. For example, if y (t) = A, sin ( 27ij\t) is phase-shifted 

by as much as n / 2, where A, and /, are the amplitude and frequency, respectively, the 
received signal is t](t) = A x sin (2^// + n! 2) = A x cos(2^/jt) which is orthogonal to 
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.s’, (t ). The integrator output is then zero and the signal cannot be detected correctly. In 

such cases, the MFSK signals can then be recovered using noncoherent detection , which 
requires only a match in frequency between transmitted and received signals. Two 
noncoherent demodulation schemes, matched filter envelope detector and quadrature 
correlator square law detector are described here [6]. 

The matched filter envelope detector consists of two stages as shown in Figure 2. 
The channel has phase-shifted the transmitted signal s(t) by 6 radians. The received 
signal s(l,6) + n(t ), where n(t) is the channel noise, enters the matched filter at the first 
stage. The matched filter has an impulse response that is the mirror image (time-reversed 
version) of the transmitted signal s(t ) delayed by the symbol duration T s . This 

maximizes its output signal-to-noise ratio when the frequency f\ of the carrier signal is 

matched. At the second stage, the envelope detector consisting of a rectifier and low pass 
filter is then used to extract the signal's complex envelope from the matched filter output. 
At the decision stage, the maximum detector decides the symbol transmitted S t 

(i = l,...,M) by selecting the largest decision statistic X i at each sampling period. 

Finally, the symbol is mapped into bits to recover the digital information. Note that the 
phase information of the signal was not used to define the matched filter or the envelope 
detector. 



Figure 2. Block diagram of noncoherent MFSK demodulation using matched filter 

envelope detection. 
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The quadrature correlator square law detector shown in Figure 3 uses both the in- 
phase ( 1 ) and quadrature ( Q ) channels for detecting the carrier signal frequency f] in 

each of the M detection branches. The setup is similar to a correlator in the coherent 
detector, but mitigates the effect of phase shifts. Squaring the correlator output removes 
the negative values and gives the magnitude. This has a similar function to the envelope 
detector. For example, if signal s x (t) was transmitted and there was no phase shift 
(6* = 0 radians), the maximum correlator output will appear at the in-phase channel at the 
first detection branch z u 2 ■ If .v, (t) has phase-shifted by tt / 2 radians, the maximum 

output appears at z x Q 2 . The remaining detection branches should give near zero outputs 

for orthogonal signals such as MFSK. This maximizes the decision statistic X 2 for s x (t) 
to be detected correctly even with a phase shift of 6 radians. 



x M (f,J r/2) 


Figure 3. Block diagram of noncoherent MFSK demodulation using quadrature 

correlator square law detectors. 
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The BER for a MFSK system using noncoherent detection in AWGN from Table 
7.9 in [5] is 


[ b , noncoheren t MFSK 


M 

= —e 
4 


log 2 M 


(7) 


where M is the number of symbols, E b / N a is the bit energy-to-noise ratio and is 
related to the signal-to-noise ratio per symbol SNR and symbol energy E s by 


SNR = ^-= M . (8) 

The BER performance of noncoherent MFSK and CSK in AWGN using Equation 
(7) for M =2, 4, 8, 16, 32 and 64 is shown in Figure 4. As E b / N a increases, the BER 
decreases correspondingly, resulting in waterfall-shaped graphs. This indicates that the 
performance of MFSK and CSK improves at higher E h / N a [7]. Additionally, as the 
order of modulation M increases, the performance of MFSK and CSK improves with 
higher M , giving curves with sharper drops in BER at higher values of E h / N a in 

logarithmic scale. As a consequence of this improvement, however, the bandwidth 
required by the system must also be increased. 
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Figure 4. Bit error probability for MFSK and CSK in AWGN channel (From [8]). 


A. BIT ERROR PROBABILITY OF MFSK IN RAYLEIGH FADING 
CHANNEL 

Considering the effect of channel tap \h\ of a Rayleigh fading channel on the bit 

error probability of noncoherently detected MFSK, we obtain the conditional bit error 
probability for noncoherent demodulation of MFSK as [5] 


2 M/2 

~~M -- 


1 k=l 


M -1 

V * J 


1 e -{i^F sm 


k +1 


(9) 


where M is the modulation order. The channel tap |/i| has a Rayleigh distribution as 

f\,Xx) = 2xe~’ r , (10) 
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and the channel tap h is normalized such that E 



= 1 . 


The closed form expression 


of the corresponding MFSK’s BER for a Rayleigh fading channel from Equation (9) can 
be expressed as 


Pu = 


M/2 


M- 


-) M -1 

I 

1 k=\ 


M -1 


l + k(l + SNR) ' 


( 11 ) 


In the following section, the effect of Doppler shift on the performance of MFSK 
in Rayleigh fading channel will be considered. 


B. DOPPLER TRACKING OF MFSK 


To study the effects of Doppler shift on noncoherent demodulation, consider a 
MFSK signal with a tone of f k =k/T s . It has a Doppler shift of f d and fade h k from the 
propagating through the slow fading channel and is represented at the receiver as 
s(l) = |/z /: | /t cos[2tz-( /) ( + f d ]l + #,_] + n(t). The demodulator employs the quadrature 

correlator square-law detector, and there are l different tones in the MFSK scheme. At 
the l -th tone’s in-phase channel, the decision sample was obtained by multiplying with 
the signal ^2/T s cos (2 7rf,t ), and the result was integrated from 1 = iT s to 1 = (/ +1)T . 
The l -th tone’s /-decision sample is [5] 


x u = C' ,T ' (\K I Acos[2^(/ t + f d )t + 0 k ])yl2/T s cos 2 nf x t dt + N u 


\h,\A r(i+l)T s 




JJ 1 COS (in [(fc - /) / T s + f d ] t + e k ) dt + N u 


( 12 ) 


\h„ 


[AX sin(2 n[(k -1) + f d T s ](i + l) + 0 k )~ sin(2 n [(*-/) + f d T ]i + 6 k ) 


+ N 


l,I 


2 2n[(k-l) + f d T,] 

where I h k | is the fading channel tap, A is the signal amplitude, f k is the frequency of the 


MFSK signal, f d is the Doppler shift, 6 k is the phase of the MFSK signal and N ,, is the 
noise at the in-phase channel of the / -th tone. 
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At the /-th tone’s quadrature channel, the signal s(t)was multiplied with 
yj2/T s sin(2and integrated from t = iT s to l = (/ + !)7 . The /-th tone’s In¬ 


decision sample is 

_ r o+dt; 


JJ + ’ ’ (\K I Acos [2 n(f k + f d )t + 0 k ])yj2/T s sin 2 nf x t dt + N lQ 


JiT s 

' lK ' sin ( 2 Ah - 0 / T, + /„ ] 1 + 0,) * + N, , t 




(13) 


l A . 


j AX cos {in [(fc - /) + f d T s ] {i +1) + 0 k ) ~ cos (2 n [(fc - / ) + f d T s ] i + 0 k ) 
2 2 n[{k-l) + f d T s ] 


+ N; 


l,Q 


where N, Q is the noise at the quadrature channel of the / -th tone. In both cases, the 

higher frequency terms are dropped off as they become insignificant after filtering by the 
low-pass filter before the integrator. The terms N t Q and N, Q result from the noise and 

are modeled as independent Gaussian random variables with zero mean and variance 
<r = N o / 2. The decision variable of the l -th tone is given as 

v=A 2 ,+v J 0 - a-*) 

With the given channel tap \h k |, X, has a Rician distribution, and its probability 
density function (PDF) is 


/*(*,)= 


U 2 J 


(15) 


where a is the noise variance and I o is the modified Bessel function of the first kind of 
zero order. Furthermore, the parameter m ; 2 is the sum of the squares of the means of 
X, j and X l Q and is expressed as 


m, 2 =(Xu) 2 +(X,, Q ) 2 


= h,\ E. 


sin 


■(xfj.) 


\2 ' 


(16) 


(x[(k-l) + f„Tj 

As defined in Equation (8), the symbol energy is E s = A 2 T S / 2. At the k -th 
branch of the demodulator, we see from Equation (16) that 
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m, 


2 r -|2 

= h t “ £ I sin ( 7tf d T s ) / 7tf d T s . The conditional symbol error probability of MFSK 


is 




= i-f 

Jo 


M 

Yl\o fx,( x ') dx i 


1=1 

l*k 


M 


f Xt (x k ) dx k 


1=1 

l*k 




(17) 


fx k { x k) dx k- 


The higher frequency terms are filtered out by the low-pass filter in the demodulator and 
do not affect the decision statistics. Note that the integration in Equation (17) can be 
expressed as the Marcum Q-function, which is defined as 

Q(a,b) = f xe ■' “ ^1 (ax) dx = 1 - [ xe ^I (ax)dx 

Jb JO 


= e 


'■ 2 +b 2 )n^f a^ 


2t 

n=Q J 


(18) 


= l-e 




-( c 2 +b 2 )/2 


X “ W , b > a > 0. 
—1 \& J 


The conditional symbol error probability is now 


/j(w./,H-rn 


1=1 

l^k 


« J m, x k ^ 

1 -Q 

Vo- cry 


fx t M dx k- 


(19) 


Finally, the bit error probability for MFSK is calculated as 


n(/„)= 


M/2 


( 1 M 


M -1 


t M 


k= i 


( 20 ) 


where M is the modulation order, and the expectation of the conditional symbol error 
probability is given as 

E [P,(\K\J,]\ = \y,(\K\’fdf u (hMK ■ (2D 

The PDF of the channel tap \h k | is that of a Rayleigh distribution and takes the form in 
Equation (10). 


16 



C. NUMERICAL INTEGRATION FOR MFSK WITH MATLAB 

The computer program in Matlab (listed in Appendix 1 to 7) is used to evaluate 
the numerical integration of Equation (20) and generate the performance charts for 
MFSK with Doppler shift in the Rayleigh fading channel. The number of symbols M 
included are M =2, 4, 8, 16, 32 and 64. Each chart contains six graphs plotted with 
increasing values of f d T s = 0, 0.10, 0.20 and 0.30 with the symbol rate kept constant at 
lk symbols/sec. The symbol energy was also kept constant at 1 joule as in the case 
where the transmitter was transmitting at constant power. It was assumed there is no 
channel coding. The case when f d T s = 0 illustrates the performance of the MFSK 
scheme when there is no Doppler shift. To put the numbers into perspective, consider the 
example of an aircraft flying at a radial velocity of 300 km/h (about 83.3 m/s) towards its 
ground station while transmitting information bits at the carrier frequency of 1 GHz. 
Using Equation (1), we find the maximum experienced Doppler shift is 
f d = lxlO 9 x83.3/3xl0 s =278 Hz. Since the symbol rate was at lk symbols/sec, the 

symbol duration is T s = 1x10 3 , and the parameter f d T s is 278x10 3 = 0.28 . For a UAV 

of medium altitude long endurance (MALE) class, the cruising speed for surveillance 
missions is typically 70 knots (about 36 m/s). With a datalink transmitting at a carrier 
frequency of 2.5 GHz , the Doppler shift is f d = 2.5 xlO 9 x36/3xl0 8 = 300 Hz , and 
f d T =300x10 3 =0.30. 

The v-axis of E b / N a is chosen for the range of 20 dB to 45 dB to include the 

data points around the BER of 10 4 , which is generally the required performance for 
effective digital communications. In this range of E h / N a , the BER also exhibits a 

steady trend decrease (improved performance) as E b / N a increases. 


17 



1. The Required ^ for the Same BER Increases as Doppler Shift 
Increases 

From Figure 5, the first general observation is that as the bit energy-to-noise ratio 
E b / N a increases, the probability of bit error decreases. This was expected since for 

higher E b / N a , the bit energy is much larger than noise power, and the signal can be 

demodulated with fewer errors. For a BER of 10 4 at M = 2 with zero Doppler shift, the 
required E b / N a for MFSK in Rayleigh fading channel was approximately 40 dB, or 
27.6 dB higher as compared to 12.4 dB in AWGN (Figure 4). For the Rayleigh fading 
channel, a much higher signal-to-noise ratio is always required to achieve the same bit 
error probability as compared to AWGN only. There was a minor difference of 
approximately 0.25 dB for Doppler shifts between f d T\ = 0 and 0.10. This indicated 

that the Doppler shifts at f d T s =0.10 and below do not affect the bit error probability 
significantly. For M - 2 in the region of E b / N 0 between 20 dB and 45 dB , when the 
Doppler shift was increased from f d T s = 0.1 to 0.2, about ldB higher bit energy is 
required to achieve the same BER of 10 4 . Similarly, when f d T s was increased from 0.2 
to 0.3, the additional required E b /N 0 is now 3dB higher for the same BER. For 
M =64, the difference in E h / N a was 0.5 dB between Doppler shifts of f d T s =0.1 and 
0.2. As the Doppler shift further increased from f d T s = 0.2 to 0.3, the E b / N a required 
for BER of 10 4 was about 1.2 dB higher. 

This shows that as the amount of Doppler shift increases, the additional SNR 
required for achieving the same BER also increases. 

2. Diminishing Improvement in BER as ^ h ^ Increases for Rayleigh 

Fading Channels 

All the performance results indicate that for the same M and Doppler shift, the 
improvement in BER for MFSK became marginal as E b / N 0 increased above 20 dB in 

the Rayleigh fading channel. Because of this, the results for MFSK do not follow the 
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waterfall-shaped performance curve found in the case of AWGN where BER improved in 
proportion to the increase in E h / N 0 . Instead, the results for Rayleigh fading channel 
resembles an inverse-linear function shaped curve. The BER declines steadily (linearly 
in the log scale) even for high values of E b / N a , e.g., 45 dB . This steady trend of BER 

improvement holds true for all the orders of modulation considered from M = 2 to 64, as 
indicated in Figure 11. 

As the bit energy relative to noise power is increased, the output of the correlator 
in the matching branch to the incoming signal increases correspondingly. At the same 
time, correlator outputs from other branches are also increased in proportion due to the 
larger incoming signal. This leads to the observation that the BER is improving steadily 
as E h / N a increases. 

3. The BER Performance Improved With Diminishing Returns as Order 
of Modulation M Increased 

For the same BER, the required E b / N n is reduced with higher modulation order 

M , and this is evident from Figure 11. At the BER of 10 4 with a Doppler shift of 
f d T s =0.3, the required E b /N a decreases by about 2.5 dB from 43 dB for M = 2 to 

40.5 dB for M = 4. When the modulation is increased to M = 8, the required E h / N a is 
39.3 dB (Figure 7), which is 1.2 dB lower than M = 4. From M =8 to 16 (at 38.5 dB in 
Figure 8), the reduction in E b / N a was 0.75 dB . From M = 16 to 32 (at 38 dB in Figure 

9) , this value is reduced to 0.5 dB . Finally, from M =32 to 64 (at 37.7 dB in Figure 

10) , the reduction in E b /N a is 0.3dB. This indicates that the improvement in BER 

obtained by increasing M has a lesser effect as M gets larger. Unlike phase modulation 
schemes such as quadrature amplitude modulation (QAM) or M-ary phase shift keying 
(MPSK), there is better performance for MFSK as the order of modulation M increases. 
This improvement comes at the expense of a larger bandwidth of the system. 
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4. The E h / N a Difference With Increasing Doppler Shifts Reduced as 
Order of Modulation M Increased 


It was observed that at the same BER, the difference in E b / N a between graphs of 


different Doppler shifts is reduced as the order of modulation M increased. At M = 2 
and for a BER of 10 4 , the difference in E h / N a between graphs of Doppler shifts with 


f d T s =0.1 and 0.3 is 3 dB. At M = 4, the difference is reduced to 2.5 dB for the same 


BER. Similarly, this difference in E b / N a is further reduced from 2.2 dB for M = 8 to 


1.8 dB for M=16. With the higher order of modulation of M =32 and 64, the 
difference is almost similar at about 1.6 dB . 


This indicates that in Rayleigh fading channels with higher order of modulation 
for MFSK, the change in BER performance with Doppler shifts of up to f d T s =0.3 is 
reduced. 
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Figure 5. Bit error probability for MFSK, M = 2 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Bit error probability for MFSK, M - 4 with increasing Doppler shifts from 
f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 7. Bit error probability for MFSK, M = 8 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 8. Bit error probability for MFSK, M = 16 with increasing Doppler shifts 
from f rJ T =0 to 0.30 for a Rayleigh fading channel. 
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Figure 9. Bit error probability for MFSK, M = 32 with increasing Doppler shifts 
from f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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10. Bit error probability for MFSK, M =64 with increasing Doppler shifts 
from f d T s = 0 to 0.30 for a Rayleigh fading channel. 



20 25 30 35 40 45 


E b /N o < dB > 


Figure 11. Bit error probability for MFSK, for Doppler shift of f d T s = 0.30 with 
increasing M = 2, 4, 8, 16, 32 and 64 for a Rayleigh fading channel. 
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D. CHAPTER SUMMARY 


In this chapter, the effects of Doppler shift on noncoherently detected MFSK 
signals in a Rayleigh fading channel with AWGN were considered. As expected, larger 
E b / N a is required to offset the effects of increasing Doppler shift. In a Rayleigh fading 

channel, the BER performance does not significantly improve with higher E b / N o . 
Increasing the E h / N n only brought about gradual but diminishing improvements to the 

BER. As the order of modulation M increases, the performance of MFSK is improved 
at the cost of increased bandwidth. Additionally, the differences between the BER 
performances at increasing Doppler shifts are smaller for higher order modulation such as 
M =64. 

The same treatment and analysis will be applied to consider the BER performance 
of noncoherently detected CSK with Doppler shift in the following chapter. 
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IV. PERFORMANCE OF CODE SHIFT KEYING WITH 
NONCOHERENT DETECTION 


CSK is another form of M-ary orthogonal signaling that employs a set of 2 k 
Walsh functions to define the complex envelope of the sinusoidal signals for M distinct 
symbols. The Walsh functions are trains of square pulses obtained from the rows of the 
Hadamard matrix, which arc orthogonal to one another. The Hadamard matrix itself is 
mathematically a square matrix with all its elements as +1 or -1 and generated 
recursively as 




"1 1 1 

h m = 

HmI 2 

H Ml 2 

1 -ij 

M 

_Hm /2 

— H M 12 _ 


( 22 ) 


where M is the order of the Hadamard in powers of two. The elements in a row of the 
matrix form the chip pattern of a single Walsh function. From Equation (2.52) in [5], the 
Walsh functions w n (t) of M chips are 


M 

= n=l,2 . M (23) 

i =1 

where T is the chip duration, p(t j is a pulse with duration of 7 and unit 

amplitude, shifted by i-th chip. The term h ni refers to the (n,i) element of the matrix 


H M and has values of +1 or -1. With wjt) as the complex envelope of the signal for 

CSK with M distinct symbols, each sinusoidal Walsh signal at the baseband is 
represented from Equation (6.138) in [5] as 

s„(t) = Aw„(t)cos(2nf c t), n = l,2,...,M (24) 

where A is the amplitude and cos (2 nf c t) is the carrier of the signal. Quadrature 

correlator square-law detectors similar to that of MFSK are used to demodulate the 
Walsh signals noncoherently. In the following section the bit error rate (BER) 
performance of CSK for a Rayleigh fading channel with Doppler shift is discussed. 
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A. 


DOPPLER TRACKING OF CSK 


Consideration is now given to the effect of Doppler shift on noncoherent 
demodulation of CSK in a Rayleigh fading channel. The received CSK signal for the z-th 
symbol at carrier frequency f c with Walsh function w k (t) after experiencing a Doppler 

shift f D and channel fade \h k is s(t) = \h k | Aw t (f)cos[2;r(/ c + f d )t + 0 k j + n(t). 


For the in-phase channel of the demodulator, the signal is demodulated by first 
multiplying with yj 2 / T s w k (t) cos (2 ftf c t) and the result is integrated over one symbol 

duration from t = iT s to 1 = (z +1) 7). From Equation (10.149) in [5], the A-th Walsh 
signal’s I-decision sample is 

x u = j^ " 7 ' (| K | Aw k (0COS \ln(f c + f d )t + e k ])y]2/T s w k (t )cos 2 nf c t dt + N kJ 

(i+l)T s , x 

cos (2 nf d t + 0 k )dt + N kI (25) 

iT s v ' 

A% sin (2 nf d t\i +1) + 9 k .) - sin (2 nf d T s i + 0 k ) 

2 2 nf d T s + kJ ' 

For the quadrature channel, the received Walsh signal s(t) is multiplied with 
y]2IT s w k (t)sm(2/rf c t) and, as before, the result is integrated for one symbol duration 
from t = iT s to t = (i + l)T s . The A-th Walsh signal’s Q-decision sample is 




X k,Q = JJ ' (\K I Aw k (0 cos [ 2 x{fc +fd) t + 0 k ])P /T s w k (0 sin 2 xf c t dt + N kJ 

\ h M 


Sin ( 27 r fd f + 9 k ) dt + N U 


(26) 


K 


\ a 2 T s cos ( 27zf d t\i +1) + ^) - cos ( 27rf d T s i + 6 k ) | ^ 
V 2 2 nf d T s + ; 


For both Equations (24) and (25), the higher frequency terms are dropped as they are 
filtered out. The noise components N k , and N k Q are independent Gaussian variables 

with zero mean and variance cr 2 = N a / 2. The decision variable of the A-th Walsh 
function is 


26 



(27) 


** = 


y[Ki 


+ X 


k,Q ' 


With the given channel tap h k , X, has a Rician distribution, and its PDF is 




fm k x k^ 

2 

G 


(28) 


where I 0 is the modified Bessel function of the first kind and zero order. The parameter 
m k 2 is the sum of the squares of the means of X k 7 and X k Q and is expressed as 

m k 2 =( i w ) 2 +( x k , Q ) 2 


=h\ 2 e, 


sin 2 {xf d T s ) 


(29) 


(*/Z) 

where E s = \jA 2 T s / 2 is the symbol energy. For the other branches in the demodulator, 
the decision variable is denoted as X t for the Z-th branch that has w,it) as its Walsh 
function. When there is no Doppler shift, X , is the envelope of the Gaussian noise and is 
Rayleigh distributed. For the in-phase channel of Z-th branch, the Z-decision sample is 
given as 

x u = J^" ,7; (K | Aw k (0 COS [2 n(f c + f d )t + e k ])yj2 7¥ s w, (t) cos 2 nf c t dt + N u 
~ w* it) w, it) cos (2 7i f d t + 0 k ) dt + N u (30) 

\h, A r(i+i)r, , . 

* J j= [ T Wj it) cos ( 2jif d t + 0 k ) dt + N LI . 

In order to evaluate the integral, consider the case of the Walsh function w Xt) 


with M 12 number of +1 followed by M 12 number of -1, which is the (Ml 2 + l) ,,! 

row of the Hadamard matrix. There is no loss of generality since the product of two 
Walsh functions in the same set, i.e., w k (t) and w,(t) in this case, is another Walsh 
function. Equation (29) becomes 
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p(i'+l/2)7^ / v f(l+L)T s , V 

J. r cos (2, 7 /, 7 + 0 k ) * - J cos( 2 ^// + 6 > t ) * 


•J^T, 


a+m 

(i+l/2)T s 
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sin (2 nf d T s (i +1 / 2) + 0 A ) - sin (2 nf d T s i + 0 k ) 

2 *f d T, 

sin (2 nf d T s (i +1) + 0 k ) - sin (2 nf d T s (i + U2) + 0 k ) 
2nf d T s 

2 sin (2 nf d T s (i +1 / 2) + 9 k ) - sin (2 nf d T s i + 0 k ) 


+ N U (31) 


2 nf c ? s 

sm(2xf d T s (i + l) + 0 k ) 


+ N, 


l,I * 


2 nf d T t 

Similarly, for the /-th Q-decision sample, the signal is multiplied with 
^2/T s w,(0sin(2;r/ c t) anc j j s gj ven as f 0 n 0W s from Equation (10.156) in [5]: 

x i, Q = C"" (K I Aw k W cos [ 2 n(f c +fd)t + 0 k ])V 2 !T S w, (t) sin 2 nf c t dt + N lQ 


j= £r' ,,; Wj (0sin (2 nf d t + 0 t ) dt + IV,■ fi 


(32) 




2cos(2^/ d r s (j + l/2) + ^)-cos(2;r f d T s i + 6 k ) 
cos (27rf d T s (i +1) + 0 k ) 


2 nf d T s 

As a result, the sum of the squares of the means is 

\2 /- \2 


+ ty. 


m{ = (X/./)“+(X/. G '" 


2 r 16sin 2 (nfj's /2)-4sin 2 (xf d T s } 

k I _ „ „ \ 2 


(33) 


( 2 ^/Z) 

The BER of CSK is identical to MFSK as in Equation (20). 
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B. NUMERICAL INTEGRATION WITH MATLAB 

As in the case for MFSK, a similar Matlab program (listed in Appendix 8 to 11) 
was used to compute the BER for CSK. The main difference between the program for 
MFSK and CSK is in the function for computing the mean of the sum of squares. This 
was different for the /-th and the k- th Walsh signal’s I-decision sample, which was m 2 

and m k 2 , respectively, as defined in Equations (29) and (33). As before, M denotes the 
order of modulation and includes the values of 2, 4, 8, 16, 32 and 64. Each of the charts 
is plotted for the increasing values of f d T s = 0, 0.10, 0.20 and 0.30 with the symbol rate 
kept constant at 1 k symbols/sec. The symbol energy is kept constant at 1 joule. 

The following results and observations are noted from the charts in Figures 12 to 
19 from the numerical integration. 

1. For the Same BER, CSK Requires Higher E b / N a Than MFSK When 
There Is Doppler Shift 

Considering M = 2 with f d T s =0.3 and a BER of 10 4 , we see that the required 
E b / N a for CSK is 45 dB (Figure 12). This is approximately 2 dB higher than MFSK. 
At M = 4 the E b / N a for CSK is 43.1 dB (Figure 13) while MFSK is at 40.5 dB, and 
their difference is increased to 2.6 dB . At M = 8 CSK requires 42.4 dB (Figure 14), 
which is 3.1dB higher than the required E h /N a of 39.3 dB for MFSK. At M =16 the 
E b / N a for CSK is 42.3 dB (Figure 15) while MFSK requires 38.5 dB, which is 3.8 dB 
lower. At M =32 the E b / N a for CSK is 42.1 dB (Figure 16), while MFSK is at 38 dB 
and the difference is 4.1dB. At M = 64 the E b / N a for CSK is also approximately 
42.1 dB (Figure 17), while MFSK is at 37.7 dB and the difference is further increased to 
4.4 dB. 

It is noted that CSK consistently requires a higher E b / N o to achieve the same 
BER as MFSK for a Rayleigh channel when there is Doppler shift. For higher 
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modulation orders, the additional required E b / N a is greater, up to 4.4 dB for M =64 . 

This occurs because CSK at higher modulation orders of M uses Walsh functions with 
chip patterns that have more transitions per symbol time. For example, for a code with 
alternating +1 and -1 square pulses, i.e., [+1,-1, +1,-1... + 1,-1], there are 64 transitions 

for M = 2 b as compared to 32 transitions for M = 2'' which is halved the number. In the 
frequency spectrum, the code with more chip transitions is more susceptible to Doppler 
shift since it is completely unrecognizable by the correlator as compared to codes with 
fewer transitions. In the case of MFSK, each symbol was represented by a tone f k 

having a separation of 1 IT s from adjacent tones to ensure orthogonality. This tone 

separation is maintained even as M increases; therefore, in the frequency spectrum, 
Doppler shift affects MFSK signals to a lesser extent than CSK signals at higher 
modulation orders. Comparing the BER graphs of MFSK and CSK in Figure 19, we 
observe that MFSK at M = 32 performs better than CSK at M =64 by about 4.1 dB . In 
general, MFSK exhibited better BER performance with Doppler shift in Rayleigh fading 
channel as compared to CSK. 

2. For the Same BER, the Required E b l N a Increases as Doppler Shift 
Increases 

Similar to MFSK, noncoherently detected CSK exhibits the trend that the required 
E b / N a to achieve the same BER is larger as the Doppler shift increased. For M = 2 at 

Doppler shift of f d T s = 0.3, the CSK modulation at BER of 10 4 requires E b / N a of 
about 45 dB, which is 5 dB higher than zero Doppler shift at 40 dB. A similar 
observation is made for M = 64 at the same BER of 10 4 , which requires 42.1 dB at 
Doppler shift of f T = 0.3, or approximately 5.4 dB higher than the zero Doppler shift 
case at 36.7 dB . This shows that for CSK, the increase in E b / N 0 from a Doppler shift 
of f d T =0.1 to 0.3 is fairly consistent at 5 dB. This is different from the case of MFSK, 
where the difference in E b / N a from M = 2 to 64 decreases from 3 dB to 1.6 dB . 
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As in the case for MFSK, increasing E b / N a does not result in the BER declining 
rapidly for CSK. The BER performance of CSK gradually improved with higher E b / N n 
in an inverse-linear fashion. 

3. The BER Performance Improves With Diminishing Returns as the 
Order of Modulation Increases 

Similar to MFSK analyzed in the earlier sections, CSK is an orthogonal signaling 
scheme, and its BER performance improves as the modulation order M increases 
(Figure 18). This improvement, however, shows diminishing returns. For a BER of 10 4 
at f tl T = 0.3 and going from M =2 (Figure 12) to M =4 (Figure 13), the required 
E b /N a is reduced by the greatest amount, approximately 1.9 dB from 45 dB to 
43.2 dB. From M = 4 to 8 (E h / N o of 42.4 dB for BER of 10 4 in Figure 14) the 
reduction is 0.7 dB. By further increasing M = 8 to 16 (at 42.3 dB in Figure 15), the 
reduction of 0.1 dB is insignificant. The performance for M = 32 and 64 is similar at 
E b / N o of 42.1 dB (Figures 16 and 17), or approximately 0.2 dB lower than the case of 
AT = 16. 

This indicates that by increasing M, the improvement in BER has a lesser effect 
when M is greater than 8. At higher modulation orders, the required bandwidth of the 
system also increases since longer Walsh codes with higher number of chip transitions 
are used. 
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Figure 12. Bit error probability for CSK, M = 2 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 13. Bit error probability for CSK, M = 4 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 14. Bit error probability for CSK, M = 8 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 15. Bit error probability for CSK, M = 16 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 16. Bit error probability for CSK, M = 32 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 17. Bit error probability for CSK, M = 64 with increasing Doppler shifts from 

f d T s = 0 to 0.30 for a Rayleigh fading channel. 
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Figure 18. Bit error probability for CSK, for Doppler shift of f d T s = 0.30 with 
increasing M = 2, 4, 8, 16, 32 and 64 for a Rayleigh fading channel. 
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Figure 19. Bit error probability for MFSK and CSK, for Doppler shift of f d T s = 0.3 

with M = 32 and 64 for a Rayleigh fading channel. 
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C. CHAPTER SUMMARY 

In this chapter, the performance of noncoherent CSK was evaluated for increasing 
Doppler shift for a Rayleigh fading channel with AWGN. It was observed that for the 
same modulation order M and BER, CSK requires a higher E h / N a for the same BER 

than MFSK. With higher modulation, the number of chip transitions in the Walsh codes 
is increased. In the frequency spectrum, Walsh functions with more transitions occupy a 
much larger bandwidth than those with a lower number of transitions. Hence, the Walsh 
codes are affected by Doppler shift differently and are more susceptible, resulting in 
higher BER (worse performance). As an orthogonal signaling scheme, CSK is similar to 
MFSK in that its BER performance improves as modulation order is increased at the cost 
of larger bandwidth requirements. When the modulation order is increased beyond 8- 
CSK, this improvement in BER shows diminishing returns. 
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V. PERFORMANCE OF M-ARY FREQUENCY SHIFT KEYING 
AND CODE SHIFT KEYING IN AWGN WITH DOPPLER SHIFT 


Consider the case of noncoherently detected MFSK and CSK in the presence of 
AWGN with Doppler shift. This models the scenario of satellite communication where 
digital information is relayed between an orbiting satellite and a flying aircraft (Figure 1). 
It is assumed that there is only a dominant line-of-sight (LOS) from space to sky with no 
multipaths. Since there is relative motion between the satellite and aircraft, the 
transmitted signals are subjected to Doppler shift. 


The single channel tap |/z| of the AWGN is now set to unity, which implies that 
there is only one LOS and no attenuation due to multipath fading. The expectation term 
in Equation (20) can be expressed as P e (f d ), which is the symbol error probability as a 
function of Doppler shift alone. Equation (20) is now modified as 

Ml 2( 


p df,i) = 


1 M 


M -1 

NUMERICAL INTEGRATION FOR MFSK AND CSK WITH MATLAB 


(34) 


The Matlab program (listed in Appendix 12) which implements a similar 
algorithm for the Rayleigh fading channel earlier was used to calculate the bit error rate 
(BER) of MFSK and CSK in AWGN with Doppler shifts. The function to calculate BER 
is modified to evaluate Equation (34). The charts are plotted for both MFSK and CSK 
using bit energy-to-noise ratio E b / N a from 0 dB to 20 dB , showing the BER of up to 

10- 5 (Figures 20 to 34). The following observations in AWGN are noted and compared 
to results obtained earlier in the Rayleigh fading channel. 

1. For the Same BER, the Required E b l N a in AWGN Is Smaller Than 
for a Rayleigh Fading Channel 

It was found that the required E b / N a in AWGN to achieve the same BER was 
lower than in Rayleigh fading. Significantly higher values of SNR were required to 
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overcome the effects of Rayleigh fading. This observation holds true for both MFSK and 
CSK. For example, in AWGN for CSK with M = 2 , E b t N 0 of 19 dB is required to 

achieve a BER of 10 -4 at f.T = 0.3, which is 26 dB lower than E, / N of 45 dB for a 
Rayleigh fading channel. At a higher modulation order of M = 64, the required E b / N a 
is 29.2 dB lower in AWGN (12.8 dB ) compared to Rayleigh channel (42.1 dB). In the 
case for MFSK with M =2, E b / N n of 17.2 dB is required to achieve a BER of 10~ 4 at 
f d T s = 0.3 in AWGN, which is 25.8 dB lower than 43 dB for a Rayleigh fading channel. 
At M =64, this difference becomes 28.3 dB between E b l N n of 9.4 dB in AWGN and 
37.7 dB for a Rayleigh fading channel. 

The familiar waterfall shaped charts for MFSK and CSK performance in AWGN 
show that BER performance can be significantly improved with SNR. 

2. The BER of MFSK is Better Than CSK in AWGN With Doppler 
Shifts 

A general observation from performance for Rayleigh channels that still holds is 
that MFSK has better BER performance than CSK in AWGN. Comparing BER 
performance at 10 4 for f d T s = 0.3 in AWGN, we see that 2-CSK requires an E b / N 0 of 
19 dB while 2-FSK requires 17.2 dB , approximately 1.8 dB . Similarly, 64-CSK requires 
an E b /N a of 12.8 dB while 64-FSK requires 9.4 dB, which is 3.4 dB lower. As 

previously discussed, once the coded signals in CSK are Doppler-shifted, they become 
unrecognizable by the correlator in the demodulator, and this is especially the case for 
Walsh codes with more transitions. 

3. BER Improves With Increasing Modulation Order With Diminishing 
Returns 

The consistent trend for orthogonal signaling is that as modulation order M 
increases, the BER performance improves with diminishing returns. For example, at a 
BER of 10 4 , 2-FSK with /'('/' =0.3 requires E b ! N a of 17.2 dB while 4-FSK requires 
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14.2 dB and is 3 dB lower. For 32-FSK and 64-FSK, the required E b / N 0 are 10.2 dB 
and 9.4 dB , respectively. The required SNR for 64-FSK is only approximately 0.8 dB 
lower than 32-FSK. Similarly, for 2-CSK the E b l N a required is 19 dB compared to 

16.3 dB for 4-FSK. At 64-CSK the E b /N 0 required is reduced tol2.8dB, 
approximately 6.2 dB lower than 2-FSK. 

4. For Higher Modulation Orders, the Difference in E b / N a for Doppler 
Shifts Is Reduced for MFSK but Remains Constant for CSK 

As M increases, the difference in E h / N a between Doppler shift of f d T s =0.1 
and 0.3 for MFSK decreases from 4.3 dB (2-FSK) to 3.3 dB (64-FSK). In the case of 
CSK, this difference in E b / N a remains fairly constant at approximately 5.8 dB from 2- 

CSK to 64-CSK. This was attributed to MFSK performing better than CSK in the 
presence of Doppler shift. In other words, increasing the modulation order for MFSK 
helps to reduce the BER at higher Doppler shifts but for CSK has negligible effects. 
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Figure 20. Bit error probability for MFSK for increasing Doppler shifts of f cl T t = 0, 

0.1, 0.2 and 0.3 with M = 2 in AWGN. 



Figure 21. Bit error probability for MFSK for increasing Doppler shifts of f d T s = 0, 

0.1, 0.2 and 0.3 with M = 4 in AWGN. 
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Figure 22. Bit error probability for MFSK for increasing Doppler shifts of f d T s = 0, 

0.1, 0.2 and 0.3 with M = 8 in AWGN. 



Figure 23. Bit error probability for MFSK for increasing Doppler shifts of f d T s = 0, 

0.1, 0.2 and 0.3 with M = 16 in AWGN. 
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Figure 24. Bit error probability for MFSK for increasing Doppler shifts of f d T s = 0, 

0.1, 0.2 and 0.3 with M = 32 in AWGN. 



Figure 25. Bit error probability for MFSK for increasing Doppler shifts of /'/7 = 0, 

0.1, 0.2 and 0.3 with M = 64 in AWGN. 
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Figure 26. Bit error probability for CSK for increasing Doppler shifts of f d T s = 0, 0.1, 

0.2 and 0.3 with M = 2 in AWGN. 



Figure 27. Bit error probability for CSK for increasing Doppler shifts of f d T s = 0, 0.1, 

0.2 and 0.3 with M = 4 in AWGN. 
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Figure 28. Bit error probability for CSK for increasing Doppler shifts of f d T s = 0, 0.1, 

0.2 and 0.3 with M= 8 in AWGN. 



Figure 29. Bit error probability for CSK for increasing Doppler shifts of f d T s = 0, 0.1, 

0.2 and 0.3 with M = 16 in AWGN. 
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Figure 30. Bit error probability for CSK for increasing Doppler shifts of f d T s = 0, 0.1, 

0.2 and 0.3 with M = 32 in AWGN. 



Figure 31. Bit error probability for CSK for increasing Doppler shifts of f d T s = 0, 0.1, 

0.2 and 0.3 with M = 64 in AWGN. 
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Figure 32. Bit error probability for MFSK for Doppler shift of f d T s = 0.3 with 
increasing order of modulation from M = 2 to 64 AWGN. 



Figure 33. Bit error probability for CSK for Doppler shift of f d T s = 0.3 with increasing 

order of modulation from M =2 to 64 AWGN. 
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Figure 34. Bit error probability for MFSK and CSK for Doppler shifts of f d T s = 0.3 

with M = 32 and 64 in AWGN. 


B. CHAPTER SUMMARY 

In this chapter, the BER performances of noncoherent MFSK and CSK were 
evaluated for increasing Doppler shift in AWGN. The SNR required for the same BER is 
considerably lower in AWGN than for a Rayleigh fading channel by more than 25 dB . It 
is observed that for the same modulation order M and BER, CSK requires an E b / N a 
approximately 2 dB to 3 dB larger than MFSK. Consequently, the performance of 
MFSK in AWGN with Doppler shift is better than CSK. With higher modulation, the 
SNR required to achieve the same BER in AWGN is reduced for both MFSK and CSK 
but has a diminishing effect. As the modulation order increases, the difference in E b / N n 

between high and low Doppler shifts is reduced for MFSK but remains fairly constant for 
CSK. 
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VI. SIMULATION WITH MATLAB SIMULINK 


In this chapter, the simulation results obtained with Matlab Simulink model are 
discussed. The MFSK performances in both AWGN and for a Rayleigh fading channel 
with Doppler shift are simulated using two models built from the block sets in Simulink 
library. This allows results obtained earlier from numerical analysis to be compared to 
those generated by simulation. The symbol rate is set to 1 k symbols/sec and increasing 
Doppler shifts are introduced, ranging from zero Doppler to 100Hz, 200Hz, and 
300Hz. The simulation is carried out for modulation orders from M = 2 to 64. The 
results are presented in the following sections. 

A. SIMULATION OF MFSK FOR A RAYLEIGH CHANNEL WITH AWGN 

AND DOPPLER SHIFT 

The Simulink model for simulation of MFSK in the Rayleigh fading channel is 
shown in Figure 35. The Bernoulli Binary blockset continuously generates random bits 
(‘1’ and ‘0’) using a Bernoulli distribution. The number of bits used to form a symbol 
depends on the order of modulation selected. For example, one bit is required for 2-FSK 
and six bits are required for 64-FSK. The MFSK modulator and demodulator blocksets 
are used and implemented for simulation at baseband. Simulation at carrier modulated 
frequencies typically in the range of gigahertz are typically not used since a very high 
sampling rate is necessary to prevent simulation errors arising from an aliasing effect. A 
high sampling rate generates a large number of datapoints and slows the simulation 
considerably. Gray coding is used to minimize bit errors arising from symbols that were 
incorrectly demodulated. The MFSK is set to discontinuous phase so a minimum 
frequency separation of 1kHz (1 / 7 based on symbol rate of lk symbols/sec) for 

orthogonal signaling could be used. Perfect channel estimation is assumed, which 
allowed the Error Rate Calculation blockset to directly compare the transmitted and 
received symbols for determining the BER. As the BER approaches one in 10 s bits, the 
simulation must be run until a sufficient number of symbols have been transmitted at 
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higher bit energy-to-noise ratio E b / N a . In this case, the model is set to run until either 
one million symbols are transmitted or 100,000 symbol errors have been detected, 
whichever condition is met first. The SNR is specified in the E h / N a field of the AWGN 

blockset. Since the fading channel blocksets require a non-zero Doppler shift, the 
arbitrary value of 1 Hz for Doppler shift was used to simulate conditions of zero Doppler. 
The findings are summarized below. 



M-FSK 
D e m 0 d u I ato r 
Baseband 


Figure 35. Simulink model for MFSK in Rayleigh channel with AWGN and Doppler 

shift. 


1. BER Performance for MFSK for a Rayleigh Fading Channel 

In Equations (12) and (13) slow fading is assumed, so the single Rayleigh channel 
tap \h k \ is considered constant over a symbol duration. With the Simulink Rayleigh 

fading blockset, the implemented single fading channel tap | h k \ is not constant during a 

symbol duration. This is evident by observing the frequency spectrum of MFSK signals 
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for M =2 in Figure 36. The MFSK signals experience interference as Doppler shift is 
increased. For 2-FSK the two distinct peaks (see Figure 36a) indicating the presence of 
two frequencies at close to zero Doppler have now spread into each other at a Doppler of 
300 Hz (see Figure 36d). This is because the time coherence of a Rayleigh fading 

channel is short, and the channel tap | h k | is not time-invariant for the symbol duration. 

As a result, the MFSK signals are no longer orthogonal when there is Doppler shift. It is 
observed that the simulation predicted worse performance than the numerical analysis. In 
general, the simulation results still show the BER performance improves as SNR 
increases for E h / N a up to approximately 14 dB . These are trends consistent with those 
obtained by numerical analysis. It is also observed that as SNR increases (for 
E b / N a > 14 dB) the BER performance does not improve. The charts instead exhibited a 
tapering effect. This was due to the interchannel interference, which limited the BER 
performance. The charts for 2-FSK and 64-FSK are chosen to illustrate this discussion 
(Figures 37 and 38). 


51 



Magnitude-squared. dB Magnitude-squared, dB 



(c) 


xIO -4 



Frame: 943 Frequency (kHz) 




Figure 36. The Frequency spectrum of 2-FSK in Rayleigh fading channel with 
increasing Doppler shift of (a) 1 Hz, (b) 100 Hz, (c) 200 Hz , and (d) 

300 Hz. 
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Figure 37. Simulation results for 2-FSK in Rayleigh fading channel with increasing 
Doppler shift of (a) 1 Hz, (b) 100 Hz, (c) 200 Hz , and (d) 300 Hz. 
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Figure 38. Simulation results for 64-FSK in Rayleigh fading channel with increasing 
Doppler shift of (a) 1 Hz, (b) 100 Hz, (c) 200 Hz , and (d) 300 Hz. 
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B. SIMULATION OF MFSK IN AWGN WITH DOPPLER SHIFT 


In Chapter IV, the numerical analysis for MFSK and CSK in AWGN channel 
with Doppler shift was carried out. The Simulink model built for MFSK simulation in 
AWGN channel (in Figure 39) uses a Rician Fading Channel blockset so that the Doppler 
shift can be specified. The K-factor of this blockset, which specifies the ratio of power in 
the line-of-sight component (LOS) to the multipath components, is set to an arbitrarily 
large value of 10 10 . This allows the Rician fading to approximate a single direct LOS 
case. 



M-FSK 

Demodulator 

Baseband 


Figure 39. Simulink model for MFSK in AWGN with Doppler shift. 


1. BER Performance for MFSK in AWGN 

The results of the simulation are close to the earlier BER performance charts 
predicted by the numerical analysis. At BER of 10 4 for modulation order of M =2 with 
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Doppler shift of f d T s = 0.3, the difference in E b / N 0 is only about 0.05 dB between the 

results from simulation (Figure 20) and numerical analysis (Figure 40). With M = 64 for 
the same conditions, the difference between simulation and numerical analysis is 
approximately 0.9 dB . This shows that results from both the numerical analysis and the 
simulation agree with each other in the case for AWGN with Doppler shift. 



0 2 4 6 8 10 12 14 16 18 20 


E b /N o (dB) 

Figure 40. Simulation bit error probability for MFSK M = 2 with Doppler shifts from 

1 Hz to 300 Hz in AWGN. 
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Figure 41. Simulation bit error probability for MFSK M =4 with Doppler shifts from 

1 Hz to 300 Hz in AWGN. 



Figure 42. Simulation bit error probability for MFSK M =8 with Doppler shifts from 

1 Hz to 300 Hz in AWGN. 
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Figure 43. Simulation bit error probability for MFSK M = 16 with Doppler shifts from 

1 Hz to 300 Hz in AWGN. 



Figure 44. Simulation bit error probability for MFSK M = 32 with Doppler shifts from 

1 Hz to 300 Hz in AWGN. 


57 


















































































































DC 

LU 

CD 


_Q 

03 

-Q 

O 

CL 

O 



E b /N o < dB > 


Figure 45. Simulation bit error probability for MFSK M = 64 with Doppler shifts from 

1 Hz to 300 Hz in AWGN. 
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Figure 46. Simulation bit error probability for MFSK M =2, 4, 8, 16, 32 and 64 with 

Doppler shifts of 300 Hz in AWGN. 
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C. CHAPTER SUMMARY 

In this chapter, the simulation results for BER performances of noncoherent 
MFSK with Doppler shift were discussed. Two Simulink models were used: the first 
model illustrates MFSK performance in Rayleigh fading with AWGN, while the second 
is for MFSK in AWGN channel. It was found that the Rayleigh fading channel blockset 
was not time-invariant within a symbol duration, resulting in two observations. First, the 
BER performance predicted by the simulation is worse than the result from the numerical 
analysis. Second, a tapering-off effect is observed for E h /N n greater than 14 dB, 
indicating that the BER performance is limited by interchannel interference. In general, 
for E b / N a less than 14 dB, the trends of the charts from the simulation agree with the 

numerical analysis, which shows that as SNR increases, BER performance improves in a 
inverse-linear fashion. For the case of AWGN, the simulation results are close to the 
results obtained by numerical analysis. 
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VII. THESIS SUMMARY AND CONCLUSION 


The Matlab program for the numerical analysis of noncoherent demodulation of 
MFSK and CSK in both Rayleigh fading channel and AWGN with Doppler shift was 
developed in this thesis. In addition, the Simulink models for MFSK for both Rayleigh 
fading channels and AWGN were presented. It was found that the channel model of the 
Rayleigh fading blockset was time-varying within a symbol duration. This resulted in a 
worse than predicted performance by the simulation when compared to numerical 
analysis. Furthermore, BER performance tapered off at E h / N a > 14 dB due to 
interchannel interference. The Simulink model for AWGN gave results closer to the 
numerical analysis. 

Depending on the order of modulation, we found that MFSK generally performed 
better than CSK. This difference was approximately 2 dB to 4.4 dB at a BER of 10 4 
for a Rayleigh fading channel and about 1.8 dB to 3.4 dB in AWGN when considering a 
Doppler shift of f d T s =0.3. This occurs because Walsh coded signals in CSK become 
unrecognizable at the demodulator when they experience Doppler shift. This was 
especially the case for the higher modulation orders like 64-CSK, which has signals using 
Walsh codes with a higher number of transitions. 

It was also possible to improve BER performance for noncoherent demodulation 
of MFSK and CSK by increasing the order of modulation. However, it is not 
recommended to go higher than M = 8 since the improvement is marginal due to 
diminishing returns as M increases. 

A much higher SNR is required for a Rayleigh fading channel to achieve the same 
BER in AWGN. The case was presented for a BER of 10 4 and a Doppler shift 
f d T s =0.3. For CSK the difference in SNR is approximately 26 dB to 39 dB higher for 
a Rayleigh fading channel. Similarly, MFSK required a SNR of nearly 17 dB to 28 dB 
higher for a Rayleigh fading channel. 
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To improve BER performance when large Doppler shifts are expected, MFSK is 
recommended over CSK for noncoherent demodulation if bandwidth of a lesser concern 
in the mobile computing and Rayleigh fading environment considered. 
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APPENDICES - MATLAB SOURCE CODES 


A. MFSKPLOTEBNOM64.M (MFSK) 


This is the main program for MFSK performance in Rayleigh fading channel with 
increasing Doppler shift using f d T s = 0, 0.1, 0.2 and 0.3. This program generates 


charts for bit error probability P b against E h / N 0 for M = 64 , and the code is similar for 


M = 2, 4, 8, 16, 32 and 64. 


close all 

2- kkkkkkkkkkkkkkkkkkkkkk 
o 

% MFSK in AWGN Channel 

Sr kkkkkkkkkkkkkkkkkkkkkk 
O 

clear 

clc 


M=64; 

Rs=le3; 
Ts=l/Rs; 

S- kkkkkkkkkkk 
o 

% NEW CODE 

Sr ★★★★★★★★★★★ 

O 

sigma=l; 
No=sigma A 2*2; 


% symbol rate at lk symbols/sec 
% symbol duration 


% normalise sigma to 1, then calculate for Es 


pts=6; 

EbNo_dB_lower = 20; 

EbNo_dB_upper = 45; 

EbNo_dB=linspace(EbNo_dB_lower,EbNo_dB_upper,pts) 
EbNo=10. A (EbNo_dB/10); 

EsNo=EbNo.*log2(M); 

Es = EsNo*No 


o 

% M=64, fdTs=0 (fd=0Hz) 

Sr kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
o 

fdTs=0; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf( 1 Calculating for M=%d,fdTs = %l.2f. . . 1 ,M,fdTs)); 

rayleighNumMFSK_M64fdTsO=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
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tic 

rayleighNumMFSK_M64fdTsO(j)=rayleighmfskPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toc(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 
save(rayleighMFSKsavemydata(M,fdTs,tStart), 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumMFSK_M64 fdTs 0' ); 
display( 'Calculation completed! ' ) ; 

o 

% M=64, fdTs=0.1 (fd=100Hz) 

o 

fdTs=0.1; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 

rayleighNumMFSK_M64fdTs01=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

rayleighNumMFSK_M64fdTs01(j)=rayleighmfskPb(M,Es(j) , fd, Ts, sigma, EbNo_dB(j)) ; 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 
save(rayleighMFSKsavemydata(M,fdTs,tStart), 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumMFSK_M64 fdTs 01' ); 
display( 'Calculation completed! ' ) ; 

Qr-k'k-k'k-k'k-k'k'k-k'k'k-k-k-k'k-k'k-k'k'k-k-k'k-k'k-k'k-k-k 

o 

% M=64, fdTs=0.2 (fd=200Hz) 

Qr-k'k-k'k-k'k-k-k-k'k-k'k-k'k-k-k-k'k-k-k-k'k-k-k-k'k-k'k-k'k 

o 

fdTs=0.2; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d,fdTs=%l.2f. . . ' ,M,fdTs)); 

rayleighNumMFSK_M64fdTs02=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

rayleighNumMFSK_M64fdTs02(j)=rayleighmfskPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 
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save(rayleighMFSKsavemydata(M, fdTs, tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumMFSK_M64 fdTs 02' ); 
display( 'Calculation completed!' ); 

8?kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

o 

% M=64, fdTs=0.3 (fd=300Hz) 

Qrkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

o 

fdTs=0.3; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 

rayleighNumMFSK_M64fdTs03=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

rayleighNumMFSK_M64fdTs03(j)=rayleighmfskPb(M,Es(j),fd,Ts,sigma,EbNo_dB( j ) ) ; 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d,fdTs=%l.2f...' ,M,fdTs)); 
save(rayleighMFSKsavemydata(M, fdTs, tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumMFSK_M64fdTs03' ); 
display( 'Calculation completed! ' ) ; 


figure 

semilogy(EbNo_dB,rayleighNumMFSK_M64fdTsO, 

EbNo_dB,rayleighNumMFSK_M64fdTs01, EbNo_dB,rayleighNumMFSK_M64fdTs02, 
EbNo_dB,rayleighNumMFSK_M64fdTs03, EbNo_dB,rayleighNumMFSK_M64fdTs04); 
hold on 
grid on 

title(sprintf (' BER for Noncoherent MFSK in RAYLEIGH, M=%d',M)); 
xlabel( 'Eb/No (dB)' ); 
ylabel( 'Pb' ) ; 

axis([EbNo_dB_lower EbNo_dB_upper le-6 leO]); 
load( 'theory_awgn_rayleigh.mat' ) 

semilogy(EbNo_dB_theory,BER_AWGN_theory_M64, 'o'); 

semilogy(EbNo_dB_theory,BER_Rayleigh_theory_M64, '+') ; 
legend(sprintf (' fdTs=%l.2f, fd=%d Hz ',0,0),... 
sprintf(' fdTs=%l.2f, fd=%d Hz' ,0.1,100),... 
sprintf(' fdTs=%l.2f, fd=%d Hz' , 0.2,200), . . . 
sprintf(' fdTs=%l.2f, fd=%d Hz' , 0.3, 300), . . . 

sprintf( 'AWGN theory for M=%d ', M),sprintf (' RAYLEIGH theory for 
M=%d' ,M)); 
hold off 
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B. RAYLEIGHMFSKPB.M (MFSK) 


This is the function definition for P b in Equation (20) for Rayleigh fading 

channel. It was found that integrating from zero to positive infinity generated errors 
when numerical values becomes small and approaches zero. Therefore the program 
evaluates the upper limit for the integration such that the function P e still has a value of 

not lower than 10 5 . To do so, it steps through the range of values using a predetermined 
stepsize which changes at different range of E b / N a values. This allowed the limits to be 
evaluated in a shorter amount of time. 


function y = Pb(M, Es, fd, Ts, sigma,EbNo_dB) 
mySum=0; 

for k=l:M % summation 

h = @(hk)((rayleighmfskPe(k,M,Es,Ts,fd,sigma,hk,EbNo_dB)) .* 

raylpdf(hk, sqrt(.5))); 

if EbNo_dB <11 

stepsize=0.1; 

%disp('0-10dB'); 

else if EbNo_dB <21 

stepsize=0.05; 

%disp('10-20dB'); 

else if EbNo_dB <31 

stepsize=0.025; 

else if EbNo_dB <41 

stepsize=0.0125; 

else if EbNo_dB <51 

stepsize=0.005; 

else if EbNo_dB <61 

stepsize=0.001; 

end 

end 

end 

end 

end 

end 
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upper=stepsize; 

prev_upper=0; 
current_upper=h(upper) ; 

% move into the threshold first 
while current_upper<le-5 

upper=upper+stepsize; % keep advancing to the right 
prev_upper=current_upper; 
current_upper=h(upper) ; 

end 

% now move until current becomes smaller than prev 
while current_upper > prev_upper & current_upper>le-5 
upper=upper+stepsize; % keep advancing to the right 
current_upper=h(upper) ; 

if isnan(current_upper) % if value becomes NaN 

upper=upper-(2*stepsize); % back track 

current_upper=le-6; % set a value smaller than le-5 to 

exit 

end 

end 

Qx = quad(h,0,upper); 
mySum=mySum + Qx; 

end 

y=(M/2)/(M-l) * (1/M) * mySum; 


end 


C. RAYLEIGHMFSKPE.M (MFSK) 


This is the function definition for conditional symbol error probability P e in 
Equation (19) for MFSK in Rayleigh fading channel. Note that functions handlers h and 
J\ perform different tasks here. Instead of integrating to positive infinity, h is used to 
evaluate the required upper limit for the integration such that the integrand value is not 
lower than 1x10 x0 which may result in quad.m returning errors for values close to zero. 
The functions h { and Q x (j) are defined in such a manner as to avoid numerical 
limitations when Matlab runs out of bits to represent the mantissa of very small values. 
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function Qx=Pe(k,M,Es,Ts,fd,sigma,hk,EbNo_dB) 

n=length(hk); 

Qx=zeros(size(hk)); 

for j=l:n 

mk = Ml(hk(j),Es,fd,Ts,k,k); 

h = @ (xk) (prodlntQ (xk, hk ( j ), Es, fd, Ts, k, M, sigma) .* fxk (xk, mk, 
sigma)); 

hi = @(xk) ( (prodlntQ(xk,hk(j),Es,fd,Ts,k,M,sigma)-1) .* fxk (xk, 

mk, sigma)); 

stepsize = mk/10; % mk is the center, use it to define 

stepsize 

upper = mk+stepsize; 
lower=mk-stepsize; 
upperValue=h(upper); 
lowerValue=h(lower); 

% To find the upper limit of integration 

while upperValue > le-80 % keep moving to the right 
upper=upper+stepsize; 
upperValue=h(upper); 

if isnan(upperValue) % if value becomes NaN 

upper=upper-(stepsize); % back track 1 step 
upperValue=le-100; % set to very small to exit 

end 

end 

while lowerValue > le-80 % keep moving to the right 
lower=max(0,lower-stepsize) ; 
lowerValue=h(lower) ; 

if isnan(lowerValue) % if value becomes NaN 

lower=lower+(stepsize); % back track 1 step 
lowerValue=le-100; % set to very small to exit 

end 

end 

Qx(j) = -(quad(hi,lower,upper)); 


end 
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D. FHK.M (MFSK AND CSK) 


This is the function definition for the PDF of the Rayleigh fading channel tap | h k | 
which follows the Rayleigh distribution from Equation (10.56) in [5]. The channel tap is 


normalized to E 


\ h f 


= 1 . 


function y=fhk(hk) 
y = 2*hk.*exp(-hk. A 2); 


E. FXK.M (MFSK AND CSK) 


This is the function definition for the PDF of the Z-th tone decision sample X ,, 
which has a Rician distribution as shown in Equation (15). The built-in Matlab function 
“besseli” is used to implement I o , the modified Bessel function of the first kind and zero 

order. When mf the sum of the squares of the means is large, the Rician PDF fails since 
I o approaches the value of unity while the exponential term becomes very small. When 
mf is large, the PDF starts to resemble the Normal distribution, and the function 
“normpdf ’ is used as an approximation to the Rician PDF for mf > 20 . 

function y=fxk(xk, mk, sigma) 

% when the mean value is large, use normalpdf as approximation to 
ricepdf in order to avoid computational errors of exp(-inf) and 
besseli() 
if mk>20 

y = normpdf(xk,mk,sigma); 

else 

y = xk./(sigma A 2) .* exp(-(xk. A 2+mk A 2)/(2*sigma A 2)) .* besseli(0, 

(mk*xk)./(sigma A 2)); 

end 
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F 


ML.M (MFSK) 


This is the function definition for m 2 , the sum of the squares of the means of X , ; 
and X, Q as defined in Equation (16) for MFSK. 


function y=Ml(hk,Es,fd, Ts, k, 1) 


if fd == 0 

if k==l % this becomes a sine function 
y=sqrt(hk A 2 * Es); 
else 
y=0; 

end 

else if (fd*Ts)==l % to take care of the case when fdTs=l 
y=0; 

else 

y = sqrt((hk A 2) * Es * (sin(pi*fd*Ts)/(pi*((k-1)+fd*Ts))) A 2); 

end 

end 

end 


G. PRODINTQ.M (MFSK AND CSK) 


This is the function definition for the term 


/=i 
i*k 


n 1 j. Tv ( v /)^ 


which is the 


integrand in Equation (17) used to evaluate the conditional symbol error probability P e . 


function y=prod!ntQ(xk,hk,Es,fd,Ts,k,M,sigma) 


n=length(xk); 

y=ones(size(xk)) ; % let y be the value of product series, now reset to 

1 

for j=l:n 

for 1=1:M 
if l~=k 

ml = Ml(hk, Es, fd, Ts, k, 1) ; 

y(j) = y(j) * (1-marcumq(ml/sigma, xk(j)/sigma)); 
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H. RAYLEIGHCSKPLOTEBNOM64.M (CSK) 


This is the main program for CSK performance in Rayleigh fading channel with 
increasing Doppler shift using f d T s = 0, 0.1, 0.2 and 0.3. This program generates 


charts for bit error probability P h against E h / N a for M = 64, and the code is similar for 


M = 2, 4, 8, 16, 32 and 64. 


close all 

2- kkkkkkkkkkkkkkkkkkkkkk 
o 

% CSK in AWGN Channel 

2- kkkkkkkkkkkkkkkkkkkkkk 
o 

clear 

clc 


M=64; 

Rs=le3; 
Ts=l/Rs; 

2- kkkkkkkkkkk 
o 

% NEW CODE 

2- kkkkkkkkkkk 
o 

sigma=l; 
No=sigma A 2*2; 


% symbol rate at lk symbols/sec 
% symbol duration 


% normalise sigma to 1, then calculate for Es 


pts=6; 

EbNo_dB_lower = 20; 

EbNo_dB_upper = 45; 

EbNo_dB=linspace(EbNo_dB_lower,EbNo_dB_upper,pts) 
EbNo=10. A (EbNo_dB/10); 

EsNo=EbNo.*log2(M) ; 

Es = EsNo*No 


o 

% M=64, fdTs=0 (fd=0Hz) 

8~kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

o 

fdTs=0; 

fd=fdTs/Ts; % Doppler frequency shift 
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display(sprintf (' Calculating for M=%d,fdTs=%l.2f...' ,M,fdTs)); 

rayleighNumCSK_M64fdTsO=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

rayleighNumCSK_M64fdTsO(j)=rayleighcskPb(M,Es(j) , fd, Ts, sigma, EbNo_dB(j)) ; 
toe 

end 

toc(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d,fdTs=%l.2f...' ,M,fdTs)); 
save(rayleighCSKsavemydata(M, fdTs, tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumCSK_M64fdTs0' ); 
display( 'Calculation completed! ' ) ; 


2 - 

o 

% M=64, fdTs=0.1 (fd=100Hz) 

Qr-k'k-k-k-k'k-k-k-k-k-k-k'k-k-k-k'k'k-k-k-k'k-k-k-k'k-k-k-k'k 

o 

fdTs=0.1; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d,fdTs=%l.2f...' ,M,fdTs)); 

rayleighNumCSK_M64fdTs01=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

rayleighNumCSK_M64fdTs01(j)=rayleighcskPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d,fdTs=%l.2f...' ,M,fdTs)); 
save(rayleighCSKsavemydata(M, fdTs, tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumCSK_M64fdTs01' ); 
display( 'Calculation completed! ' ) ; 


o 

% M=64, fdTs=0.2 (fd=200Hz) 

Sr'k'k-k'k-k-k-k-k-k'k-k'k'k'k-k-k-k'k-k-k-k'k-k-k-k'k-k-k-k'k 

o 

fdTs=0.2; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d,fdTs = %l.2f. . . ' ,M,fdTs)); 

rayleighNumCSK_M64fdTs02=zeros(size(EbNo_dB)); 
tStart=tic; 
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for j=l:pts 
tic 

rayleighNumCSK_M64fdTs02(j)=rayleighcskPb(M,Es(j),fd,Ts,sigma,EbNo_dB( j ) ) ; 
toe 

end 

toc(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs=%l.2f. . . ' , M, fdTs)); 
save(rayleighCSKsavemydata(M,fdTs,tStart), 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumCSK_M64 fdTs 02' ); 
display( 'Calculation completed! ' ) ; 


Sr-k'k-k'k-k-k-k'k'k-k-k'k'k-k-k'k'k-k'k'k'k-k'k'k-k'k-k'k'k-k 

o 

% M=64, fdTs=0.3 (fd=300Hz) 

o 

fdTs=0.3; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 

rayleighNumCSK_M64fdTs03=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

rayleighNumCSK_M64fdTs03(j)=rayleighcskPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 
save(rayleighCSKsavemydata(M,fdTs,tStart), 

'M' , 'fd' , 'EbNo_dB' , 'rayleighNumCSK_M64 fdTs 03' ); 
display( 'Calculation completed! ' ) ; 


figure 

semilogy(EbNo_dB,rayleighNumCSK_M64fdTsO, 

EbNo_dB,rayleighNumCSK_M64fdTs01, EbNo_dB,rayleighNumCSK_M64fdTs02, 
EbNo_dB,rayleighNumCSK_M64fdTs03); 
hold on 
grid on 

title(sprintf ('BER for Noncoherent CSK in RAYLEIGH, M=%d',M)); 
xlabel( 'Eb/No (dB)' ); 
ylabel( 'Pb' ); 

axis([EbNo_dB_lower EbNo_dB_upper le-6 leO]); 
load( 'theory_awgn_rayleigh.mat' ) 

semilogy(EbNo_dB_theory,BER_AWGN_theory_M64, 'o'); 

semilogy(EbNo_dB_theory,BER_Rayleigh_theory_M64, '+') ; 
legend(sprintf (' fdTs=%l.2f, fd=%d Hz ',0,0),... 
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sprintf( 'fdTs=%l.2f, 

fd=%d Hz' ,0.1,100), . . . 

sprintf( 'AWGN theory 

for M=%d',M),sprintf('RAYLEIGH theory for 

M=%d' ,M)); 


hold off 



I. RAYLEIGHCSKPB.M (CSK) 

This is the function definition for P b in Equation (20) for MFSK in Rayleigh 
fading channel. It is also applicable to CSK, but named differently here. 


function y = Pb(M, Es, fd, Ts, sigma,EbNo_dB) 
mySum=0; 

for k=l:M % summation 

h = @(hk) ((rayleighcskPe(k, M,Es,Ts,fd,sigma,hk,EbNo_dB)) . * 

raylpdf(hk, sqrt(.5))); 


if EbNo_dB <11 

stepsize=0.1; 

else if EbNo_dB <21 

stepsize=0.05; 

else if EbNo_dB <31 

stepsize=0.025; 

else if EbNo_dB <41 

stepsize=0.0125; 

else if EbNo_dB <51 

stepsize=0.005; 

else if EbNo_dB <61 

stepsize=0.001; 

end 

end 

end 

end 

end 

end 


upper=stepsize; 

prev_upper=0; 
current_upper=h(upper) ; 
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% move into the threshold first 
while current_upper<le-5 

upper=upper+stepsize; % keep advancing to the right 
prev_upper=current_upper; 
current_upper=h(upper) ; 

end 

% now move until current becomes smaller than prev 
while current_upper > prev_upper & current_upper>le-5 
upper=upper+stepsize; % keep advancing to the right 
current_upper=h(upper) ; 

if isnan(current_upper) % if value becomes NaN 

upper=upper-(2*stepsize); % back track 

current_upper=le-6; % set a value smaller than le-5 to 

exit 

end 

end 

Qx = quad(h,0,upper); 
mySum=mySum + Qx; 

end 

y=(M/2)/(M-l) * (1/M) * mySum; 


end 


J. RAYLEIGHCSKPE.M (CSK) 

This is the function definition for the conditional symbol error probability P e for 
CSK in Rayleigh fading channel. It is the same as Equation (19) for MFSK. 

function Qx=Pe(k,M,Es,Ts,fd,sigma,hk,EbNo_dB) 

n=length(hk); 

Qx=zeros(size(hk)); 

for j=l:n 

mk = Ml(hk(j),Es,fd,Ts,k,k); 

h = @ (xk) (prodlntQ (xk, hk ( j ), Es, fd, Ts, k, M, sigma) .* fxk (xk, mk, 
sigma)); 

hi = @(xk) ((prodlntQ(xk,hk(j),Es,fd,Ts,k,M,sigma)-1) .* fxk(xk, 

mk, sigma)); 
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stepsize = mk/10; % mk is the center, use it to define 

stepsize 

upper = mk+stepsize; 
lower=mk-stepsize; 
upperValue=h(upper) ; 
lowerValue=h(lower) ; 


% To find the upper limit of integration 

while upperValue > le-80 % keep moving to the right 
upper=upper+stepsize; 
upperValue=h(upper); 

if isnan(upperValue) % if value becomes NaN 

upper=upper-(stepsize); % back track 1 step 
upperValue=le-100; % set to very small to exit 

end 

end 

while lowerValue > le-80 % keep moving to the right 
lower=max(0,lower-stepsize) ; 
lowerValue=h(lower) ; 

if isnan(lowerValue) % if value becomes NaN 

lower=lower+(stepsize); % back track 1 step 
lowerValue=le-100; % set to very small to exit 

end 

end 

Qx(j) = -(quad(hi,lower,upper)); 


end 


K. ML.M (CSK) 


This is the function definition for m 2 and m 2 as defined in Equations (29) and 
(33) for CSK. 

function y=Ml(hk, Es,fd,Ts,k,1) 
if fd == 0 

if k==l % this becomes a sine function 
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y=sqrt(hk A 2 * Es); 
else 
y=0; 

end 


else if (fd*Ts)==l % to take care of the case when fdTs=l 
y=0; 


else 

if k==l 

y = sqrt((hk A 2) 

else 

y = sqrt((hk A 2) 

sine(fd*Ts) A 2)); 

end 


* Es * sine(fd*Ts) A 2); % eqn (10.153) 

* Es * (sine(fd*Ts/2) A 2 - 


end 

end 

end 


L. AWGNMFSKPLOTEBNOM64.M (MFSK AND CSK IN AWGN) 


This is the main program for MFSK performance in AWGN with increasing 
Doppler shift using f d T s = 0, 0.1, 0.2 and 0.3. The program is similar for CSK and in 

both cases it generates charts for bit error probability P b against E b / N 0 for M - 64, and 
the code is similar for M = 2, 4, 8, 16, 32 and 64. 

2- kkkkkkkkkkkkkkkkk 
o 

% AWGN Channel 

o 

clear 

clc 

M=64; 

Rs=le3; % symbol rate at lk symbols/sec 

Ts=l/Rs; % symbol duration 

2- *********** 
o 

% NEW CODE 

2- kk-kk-kk-kk-kk-k 
o 

sigma=l; % normalise sigma to 1, then calculate for Es 

No=sigma A 2*2; 

pts=21; 
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EbNo_dB_lower = 0; 

EbNo_dB_upper = 20; 

EbNo_dB=linspace(EbNo_dB_lower, EbNo_dB_upper, pts) 
EbNo=10. A (EbNo_dB/10) ; 

EsNo=EbNo.*log2(M); 

Es = EsNo*No 


Sr-k-k-k'k'k-k'k-k-k-k'k'k-k'k-k'k-k-k-k'k'k'k-k'k'k-k-k'k-k-k 

o 

% M=64, fdTs=0 (fd=0Hz) 

Q'-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

o 

fdTs=0; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d,fdTs = %l.2f. . . ' ,M,fdTs)); 

awgnNumMFSK_M64fdTsO=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

awgnNumMFSK_M64fdTsO(j)=awgnPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toc(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs = %l.2f. . . ' ,M,fdTs)); 
save(awgnsavemydata(M,fdTs, tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'awgnNumMFSK_M64fdTs0' ); 
display( 'Calculation completed! ' ) ; 


o 

% M=64, fdTs=0.1 (fd=100Hz) 

o 

fdTs=0.1; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf( 1 Calculating for M=%d,fdTs=%l.2f. . . ' ,M,fdTs)); 

awgnNumMFSK_M64fdTs01=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

awgnNumMFSK_M64fdTs01(j)=awgnPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d, fdTs=%l.2f. . . ' ,M,fdTs)); 
save(awgnsavemydata(M,fdTs, tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'awgnNumMF S K_M 64fdTs01' ); 
display( 'Calculation completed! ' ) ; 
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9-'5 | r9r'j | r9r'j | r9r9c’9r9r9c’9c’9r9c’9r9r9r9r9r9c’9r9r9r9c’9r9c’9r9c’9r'5 l r9r 

o 

% M=64, fdTs=0.2 (fd=200Hz) 

2'9r9c9r9e9f9r9r9e9f9c9f9c9r9r9r9c9f9c9r9e9f9c9r9c9f9c9r9c9f9r 

o 

fdTs=0.2; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d,fdTs = %l.2f. . . ' ,M,fdTs)); 

awgnNumMFSK_M64fdTs02=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

awgnNumMFSK_M64fdTs02(j)=awgnPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toc(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d,fdTs=%l.2f...' ,M,fdTs)); 
save(awgnsavemydata(M,fdTs,tStart) , 

'M' , 'fd' , 'EbNo_dB' , 'awgnNumMFSK_M64fdTs02' ); 
display( 'Calculation completed! ' ) ; 


2-9r9e'j | c , 9c9r9e'j | r9c9r9c'5 l r9c'jk’9e'* , 9e9r9e'* , 9c7k’9e'5 l r9e9r9e'5 l r9c9r9r 

o 

% M=64, fdTs=0.3 (fd=300Hz) 

Qr'k'k-k'k-k-k-k'k-k'k-k'k-k-k-k'k'k'k-k'k-k-k-k'k-k'k-k'k-k'k 

o 

fdTs=0.3; 

fd=fdTs/Ts; % Doppler frequency shift 

display(sprintf (' Calculating for M=%d,fdTs=%l.2f...' ,M,fdTs)); 

awgnNumMFSK_M64fdTs03=zeros(size(EbNo_dB)); 
tStart=tic; 

for j=l:pts 
tic 

awgnNumMFSK_M64fdTs03(j)=awgnPb(M,Es(j),fd,Ts,sigma,EbNo_dB(j)); 
toe 

end 

toe(tStart); % stop the execution time 

display(sprintf (' Saving data M=%d,fdTs=%l.2f...' ,M,fdTs)); 
save(awgnsavemydata(M,fdTs,tStart), 

'M' , 'fd' , 'EbNo_dB' , 'awgnNumMFSK_M64fdTs03' ); 
display( 'Calculation completed! ' ) ; 


figure 

semilogy(EbNo_dB,awgnNumMFSK_M64fdTsO , EbNo_dB,awgnNumMFSK_M64fdTs01, 
EbNo_dB,awgnNumMFSK_M64fdTs02, EbNo_dB,awgnNumMFSK_M64fdTs03, 

EbNo_dB,awgnNumMFSK_M64fdTs04); 
hold on 
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grid on 

title(sprintf (' BER for Noncoherent MFSK in AWGN, M=%d',M)); 
xlabel( 'Eb/No (dB)' ); 
ylabel( 'Pb' ); 

axis([EbNo_dB_lower EbNo_dB_upper le-10 leO]); 

%axis([0 60 le-7 leO]); 


load(' theory_awgn_rayleigh.mat' ) 

semilogy(EbNo_dB_theory,BER_AWGN_theory_M64, 'o') ; 

semilogy(EbNo_dB_theory,BER_Rayleigh_theory_M64, '+') ; 
legend(sprintf (' M=%d, fdTs=%1.2f, fd=%d Hz ',0,0),... 
sprintf (' fdTs=%l.2f, fd=%d Hz' ,0.1,100), . . . 
sprintf (' fdTs=%l.2f, fd=%d Hz' , 0.2,200), . . . 
sprintf (' fdTs=%l.2f, fd=%d Hz' , 0.3, 300), . . . 

sprintf( 'AWGN theory for M=%d ', M),sprintf (' RAYLEIGH theory for 
M=%d' ,M)); 
hold off 


M. PB.M (MFSK AND CSK IN AWGN) 

For the case in AWGN, the channel tap is \h k | = 1 . 


function y = Pb(M, Es, fd, Ts, sigma,EbNo_dB) 
mySum=0; 

for k=l:M % summation 
hk=l; 

Qx = awgnPe (k, M, Es, Ts, fd, sigma, hk, EbNo_dB) ; 

mySum=mySum + Qx; 

end 

y= (M/2) / (M-l) * (1/M) * mySum; 

end 
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N, 


PE.M (MFSK AND CSK IN AWGN) 


This function definition is similar to the case for Rayleigh fading, but the 
integrand value is now set to not lower than lx 1(T 10 , yielding better results. 


function Qx=Pe(k,M,Es,Ts,fd,sigma,hk,EbNo_dB) 

n=length(hk); 

Qx=zeros(size(hk)); 

for j=l:n 

mk = Ml(hk(j),Es,fd,Ts,k,k); 

h = @ (xk) (prodlntQ (xk, hk ( j ), Es, fd, Ts, k, M, sigma) .* fxk (xk, mk, 
sigma)); 

hi = @(xk) ( (prodlntQ(xk,hk(j),Es,fd,Ts,k,M,sigma)-1) .* fxk (xk, 

mk, sigma)); 

stepsize = mk/10; % mk is the center, use it to define 

stepsize 

upper = mk+stepsize; 
lower=mk-stepsize; 
upperValue=h(upper); 
lowerValue=h(lower); 

% To find the upper limit of integration 

while upperValue > le-10 % keep moving to the right 
upper=upper+stepsize; 
upperValue=h(upper); 

if isnan(upperValue) % if value becomes NaN 

upper=upper-(stepsize); % back track 1 step 
upperValue=le-100; % set to very small to exit 

end 

end 

while lowerValue > le-10 % keep moving to the right 
lower=max(0,lower-stepsize) ; 
lowerValue=h(lower) ; 

if isnan(lowerValue) % if value becomes NaN 

lower=lower+(stepsize); % back track 1 step 
lowerValue=le-100; % set to very small to exit 

end 

end 

Qx(j) = -(quad(hi,lower,upper)); 
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end 


O. AUTOMATEM64.M (SIMULATION OF MFSK IN RAYLEIGH/A WON) 

This m-file was written to automate the process of running the simulations for 
MFSK for both Rayleigh fading channels and AWGN. The file name is changed to load 
either Rayleigh or AWGN model in Simulink. For different orders of modulation, M is 
changed and the m-file shown here is for M = 64. 


clear 

open_system( 'rayleigh_template_7 5_19nov' ); 
display(sprintf( 'Entering initialisation. . . ' )) ; 

M=64; 

k=log2(M); 

seed_binary = 12345; 
seed_rayleigh = 12349; 
randseed_AWGN = 12351; 
maxNumErrs=le6; 
maxNumBits=5e6; 

pts=21; 

EbNo_dB_lower = 0; 

EbNo_dB_upper = 20; 

EbNo_dB=linspace(EbNo_dB_lower,EbNo_dB_upper,pts) 

Rs=1000; %lk symbols/s 
Ts=l/Rs; 

freqsep=1000; % freq sep of 1/Ts to preserve orthogonality 
nsamp=100; 

display(sprintf( 'Completed initialisation!' )); 

Sr-k'k-k-k-k'k-k-k-k'k-k-k'k-k-k-k-k'k-k-k-k'k-k-k-k'k-k-k-k'k 

o 

% M=64, fd=lHz 

O 

fd=l; 

display(sprintf (' fd = %d Hz', fd) ) ; 
for i=l:length(EbNo_dB) 
sTime=tic; 

currentEbNo=EbNo_dB(i); 

sim( 'rayleigh_template_75_l9nov' ); % run the simulation 
filename = sprintf (' rayleighM64fdlsimEbNo%d.mat currentEbNo); 
save(filename, 'ErrorVec'); 

display(sprintf (' Completed EbNo = %d', currentEbNo')); 
toc(sTime); % time the loops 

end 

Qr-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k'k-k-k-k-k-k-k-k'k-k'k-k-k-k-k 

o 


82 




% M=64, fd=100Hz 

Sric'klclclcjclclcicjclc'kicjclclciclcic'kiclciclciclcic'kic'k 

o 

fd=100; 

display(sprintf( 1 fd = %d Hz', fd)); 
for i=l:length(EbNo_dB) 
sTime=tic; 

currentEbNo=EbNo_dB(i); 

sim( 'rayleigh_template_75_19nov' ); % run the simulation 

filename = sprintf( 'rayleighM64fdl00simEbNo%d.mat' ,currentEbNo); 
save(filename, 'ErrorVec'); 

display(sprintf (' Completed EbNo = %d', currentEbNo')); 
toc(sTime); % time the loops 

end 

o 

% M=64, fd=200Hz 

Qr-k'k-k'k-k-k-k'k'k-k-k'k-k'k-k-k-k-k-k-k-k-k'k'k'k'k-k'k'k'k 

o 

fd=200; 

display(sprintf (' fd = %d Hz', fd)); 
for i=l:length(EbNo_dB) 
sTime=tic; 

currentEbNo=EbNo_dB(i); 

sim( 'rayleigh_template_75_19nov' ) ; % run the simulation 

filename = sprintf( 'rayleighM64fd200simEbNo%d.mat' ,currentEbNo); 
save(filename, 'ErrorVec'); 

display(sprintf (' Completed EbNo = %d', currentEbNo')); 
toc(sTime); % time the loops 

end 

o 

% M=64, fd=300Hz 

o 

fd=300; 

display(sprintf (' fd = %d Hz', fd)); 
for i=l:length(EbNo_dB) 
sTime=tic; 

currentEbNo=EbNo_dB(i); 

sim( 'rayleigh_template_75_l9nov' ); % run the simulation 
filename = sprintf( 'rayleighM64fd300simEbNo%d.mat' ,currentEbNo); 
save(filename, 'ErrorVec'); 

display(sprintf (' Completed EbNo = %d', currentEbNo')); 
toc(sTime); % time the loops 

end 
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